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Abstract 



An interaction — xQ ■ Q(l + Bt(1) ■ t(2)) is used in a shell model calculation 
for 10 Be. Whereas for B = the 2f state is two-fold degenerate, introducing 
a negative B causes an 'isovector' 2 + state to come down to zero energy at 
B = —0.67 and an S = 1 L = 1 triplet (J = + , 1 + , 2 + ) to come down to zero 
energy at B = —0.73. These are undesirable properties, but a large negative 
B is apparently needed to fit the energy of the isovector giant quadrupole 
resonance. 
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I. INTRODUCTION 



In this work, we wish to deal with a mismatch which occurs when the schematic 
quadrupole-quadrupole interaction, including an isosp in-dependent part, is used as a 
particle-hole interaction in R.P.A. calculations or is used as a particle-particle interaction in 
shell model calculations. Concerning the former, we have the reviews of Bayman [|] and of 
Bes and Sorensen which show that the Q ■ Q interaction plus pairing can explain the low 
lying 2 + vibrational states in even-even nuclei. These are also well described in the books of 
Bohr and Mottelson |3] and of Soloviev Q|. Bes, Broglia and Nilsson ||, Bohr and Mottelson 
and Suzuki and Rowe || also use the Q ■ Q interaction for high frequency modes and 
note that, in order to explain the large splitting between the isoscalar and isovector giant 
quadrupole resonances, one needs a strong isospin-dependent term Q ■ Qr{\) ■ r(2). Using 
the parametrization V = ~xQ • Q(l + Bf{\) ■ r(2)) (where r(l) • f(2) = 1 for T = 1 states 
and -3 for T = states), Bohr and Motteslon |§ state that B is equal to -3.6, very large 
indeed. Soloviev || quotes the formula = — 0.5(2A + 3) which equals -3.5 for A = 2, but 
in actual calculations he uses — 0.2(2A + 3). In Ref. it is noted however that if B = —3.6 
in a large space including AiV = 2 excitations, then if one truncates to a AN = space one 
should use a value which is much smaller in magnitude B = —0.6. More recent references 
include those of Hamamoto and Nazarewicz Jj]] and of Nojarov, Faessler and Dingfelder 
|. The latter authors made a critical study of the parameter B and concluded that it 



should have a smaller magnitude than was previously used. They use B = —2 in Ref. |10 
then compare this favoured value with other values e.g. B = —3.6 and B = —0.6. 

On the other hand, the Q ■ Q interaction has been used as a particle-particle interaction 



as well, especially by Elliott with his SU (3) scheme DTTJ. In the s — d shell this interaction is 
used to explain rotational behaviour in many nuclei e.g. 20 Ne, 22 Ne and 24 Mg. The model, 



as shown by Elliott [|11[] and by Harvey also helps explain deviations from the extreme 
rotational model due to the truncation effects in the shell model. However, Elliott [O uses 
a Q ■ Q interaction without an isovector term (i.e. with B = 0). One may well wonder what 
would happen to his scheme if we introduced a large r(l) • r(2) term. 

In this work we consider precisely this problem but we work in the p shell, where things 
are even simpler than in the s — d shell, and consider the case of 10 Be. We choose this 
nucleus because it is strongly deformed {(3 = 1.12 according to the tables of Raman et. al. 
|13| ) and also because it is an iV ^ Z nucleus. For such a nucleus we can have isovector 
transitions from the ground state which don't change the overall isospin. Such transitions 
are very important for our considerations. 
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II. CALCULATIONS 



We perform p shell calculations for states in 10 Be using the interaction 

V = - X Q-Q(1 + Bt(1)-t(2)) (1) 

with x — 0.36146. We study the behaviour of selected states for various negative values of 
B, the coefficient of the isovector Q ■ Q interaction. 

An attractive feature of the multipole-multipole interaction Lr ■ Lt is that di- 
rect particle-hole matrix elements of the form ((PbH^ 1 ) jt V (PaHJ 1 ) jt ) vanish un- 
less J = L and T = r (the exchange terms are usually taken to be zero). Thus 
if one adjusts the parameters of one (JT) mode pictured as an R.P.A. state ^ JT = 
Y.{Xph{PH- 1 ) jt - YpniPH- 1 )^ JT }, then the other modes (JT)' ^ (JT) are not af- 
fected. 

The expressions that we use for the particle-particle matrix elements are as follows: 
[j A (l)j B (2)] JT \ - X Q ■ Q[(l + Bt(1) ■ t(2))\ I {[jc(1)Jd(2)} jt -[jc(2)j D (1)} JT }: 

'2 3 A 



2*7+1 



-X[l + (*r,i - 3*r, )S][l + (-l) jc+jD+J+T P jc j D ] 

x(¥ A [Y L r L & c ] jA )(¥ B [Y L r L ¥ D ] jB )U(j A L J j D ; j c j B ) 

where we define our singly reduced matrix elements by the following convention: 

(L 31 M mi \j 2 m 2 )(^[0 L ^] j2 ) = L M 

In the above U is the unitary Racah coefficient. Note that the entire dependence on J, the 
total angular momentum of the two particles, is contained in the above U coefficient. 
The expression for the direct part of the particle-hole interaction is: 

((W) JT \V\ Uc3dY T ) = -x(*r,o + B8 T>1 )[A(2j A + l)(2j c + l)]h JjL 

x (& B [Y L r L ^ jA ] jB )(^ jc [Y L r L ^ jD ] jc ) 

In the above expression we have made the isospin dependence as explicit as possible. We 
see that for particle-hole states the T = shift is proportional to — x and the T = 1 to —xB. 
We can give a simplified derivation of the value of B using the quadrupole giant resonance 
data presented in Suzuki and Rowe's work 0. The isoscalar quadrupole resonance is at 
-^r MeV. The unperturbed energy of these giant resonances is 2hu = -^r MeV. Assuming 
a Tamm-Dancoff model, the ratio of isovector to isoscalar shifts is 

-X 63 - 82 
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A. The B = Limit 



We now consider shell model calculations for 10 Be using the interaction of Eq.(l). For 
B = we simply have the interaction — xQ • Q- In the p shell, as for any spin-isospin 
independent interaction, the states are classified by the orbital symmetry values [fx f 2 fa] 
as well as by the quantum numbers L, S and T ||14|| . The energies are given by the Elliott 



formula |TTJ for SU (3) with A = fx — f 2 and \x = f 2 — fs'- 

E = x[-4(A 2 + + + 3(A + //) + 3L(L + 1)] (2) 

where x = ff^X an d & is the harmonic oscillator length parameter (b 2 = The values 
of x an d X f° r 10 -Be are respectively 0.36146 and 0.1286. The ground state has quantum 
numbers [4 2 0]5 = 0L = 0T = 1; J = + . The first excited state is doubly degenerate: 
L = 2 S = [4 2 0]; J = 2^, 2 2 , and the excitation energy is 18%. We also consider the 
next excited states arising from two degenerate orbital symmetry states [3 3 0] and [4 1 1] . 
The other quantum numbers are L=1S=1T=1. The J values are therefore + , 1 + and 
2 + . In other words, for each orbital symmetry we have a triplet of states. The excitation 
energy is 30%- These states cannot be reached from the ground state by the Ml operator. 
There are several other states, one of which is the scissors mode state with quantum numbers 
L = 1S = 0(J= 1 + ). This state has an excitation energy of 66% and orbital symmetry 
[/] = [321]. Of particular interest is the fact that the T = 1 and T = 2 scissors mode states 
are degenerate in energy for the above interaction — xQ • Q- Note that the scissors mode is 
not the lowest 1 + state; the aforementioned L = 1S=1J = 1 + states lie lower. 



B. The Dependence on B 

The main thrust of the paper is in this section. Having noted in the introduction that 
a large and negative value of B is needed to fit the splitting of isovector and isoscalar giant 
quadrupole resonances, we will now study what happens to selected states in 10 Be when a 
finite negative B is introduced. The results are presented in Fig. 1. 

We first focus on the two J = 2 + states, which for B = are degenerate and lie lowest. 
This is not the case experimentally; the 2f state is at 3.368 MeV and the 2 2 state at 5.960 
MeV the splitting being largely due to a spin-orbit interaction. Both 2 + states have the 
same orbital symmetry as the J = + ground state [/] = [4 2 0]. We see from Fig.l that as 
B is made negative the degeneracy is removed with one state going rapidly down towards 
zero energy and the other rising in energy. The behaviour is not precisely linear, but the 
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linear approximation (which would hold if there were no admixtures of states of different 
symmetry) is remarkably good for negative 5. A linear fit to the behaviour for negative B 
is as follows: 

— = 18 - 26.935 5 

X 

with Ea vanishing at B — —0.668. 

The fact that the 2\ state comes down towards zero already is a signal of a very peculiar 
behaviour. We find even more peculiar behaviour if we look at the transition rates in Table 
/. For electric quadrupole excitations we define B(E2, e p , e n ) as the transition rate when 
effective charges e p and e n are used for the proton and neutron respectively. 

In Table I we list the isoscalar transition rate B(E2, 1, 1) and the isovector transition rate 
5(52, 1, —1). Note that the state \2a) comes down in energy as B becomes more negative 
and is excited from the ground state only by the isovector operator, whilst the state |2^) 
goes up in energy as B becomes more negative and is excited only by the isoscalar operator. 
This behaviour, which is shown in Fig.l, clearly goes against experiment. The lowest 2 + 
states in essentially all nuclei, although they may have some isovector part, are dominantly 
isoscalar. 

We now look at other selected states. The states which for B = have quantum numbers 
S— 1L — If — [3 3] and [/] = [4 11] (two degenerate configurations) lead to two sets of 
triplets with total angular momenta J = + , 1 + and 2 + . When a finite negative B is turned 
on, the J degeneracy is maintained but the degeneracy between the two sets of triplets is 
removed. Both sets come down in energy, and in the linear approximation we get: 

— = 30- 41.160151 

X 

with E c vanishing at B — —0.729, and 

— = 30 - 31.850151 

X 

with E D vanishing at B = —0.942. We show the behaviour of the state |C) as a function of 
5 in Fig.l. This figure shows clearly the linear collapse of this state as well as the isovector 
2 + state | A) as a function of negative B. 

Note that the J = + , 1 + , 2 + triplet \C) vanishes at a value of B very close to that 
for the J = 2+ state \A). The values are B = -0.729 and B = -0.668 respectively. Thus, 
care must be taken not to confuse the 2 + states of each configuration in this region of 5 
and beyond. There is a small region of B from -0.67 to -0.84 where the state J = 2\ is the 
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lowest state. But then the triplet J = Oq, 1q, 2q, although starting from a higher energy 
at B = 0, has a slope of larger magnitude than the 2\ state, and ultimately becomes the 
ground state for B < —0.84. 

We finally look at the states with orbital symmetry [4 2 1]. One of these states is the 



L=1S = 0T = 1 scissors mode state and is therefore of special interest |0J, |fL6| . Equally 
of interest is the other part of the scissors mode strength L = 1S = 0T = 2. The behaviour 
as a function of B is shown in Fig. 2 and Table II. For B = the T = 1 and T = 2 scissors 



are degenerate in energy ||i7|,|14 |, and there are four degenerate states in all for each T. As B 
is made negative two T = 1 states come down in energy and two come up. We give formulae 
for only the two extreme states -one going down the fastest and one going up the fastest 
(respectively): 

E 

-JL = 66 - 48.362151 

X 

— = 66 + 23.761ISI 

X 

By looking at the Ml rates at Table II, we see that the state J — lp T — 1 and 
J = Iq T = 2 form the scissors modes -they get all the isovector orbital strength. The 
state \D) which goes down in energy has no isovector orbital strength. In Fig. 2 we show the 
behaviour, as a function of B, of the T = 1 and T = 2 scissors modes. Note that whereas 
for B = the two are degenerate, for negative B the T = 2 strength comes below the T = 1 
strength, another peculiar result. 

We show now in Table II the isoscalar and isovector (scissors mode) orbital magnetic 
dipole rates. The transitions are to one state. From B = to B = —0.6 the isovector rate 
to T = 1 final states increases from 0.0890 fi 2 N to 0.230 pi 2 N i.e. an isovector Q ■ Q interaction 
with negative B causes the scissors mode strength to increase. Conversly, for positive B, 
the strength decreases with increasing B. It should be noted that for B = the strength 
is jj- = 0.0890//^. It should also be noted that in this limit the T = 2 scissors mode is 
degenerate in energy with the T = 1 mode, and the strength to T = 2 is n 2 N . The ratio is 
= §. The isoscalar orbital rate starts at zero for B = and increases with negative 
B to a finite albeit very small value i.e. for B = —0.6 B(M1) Ti sosca £ar orbital— 2.52 x 10 -3 fi 2 N . 

The energy-weighted Ml orbital strength also increases as B is varied. The combined 
T = 1 and T = 2 energy-weighted strength in the range > B > —0.6 is given by the 
approximate linear formula 

EWS(B) = EWS(B = 0)(1 - 1.7B) 
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III. BEYOND THE CROSS OVER REGION 



When B becomes less than —0.67, the state \A) which is a J = 2 + 'isovector state' 
becomes the ground state. It so remains in the range —0.67 > B > —0.84. For B < —0.84 
the ground state becomes a triplet 5 = 1 L = 1 J = + , 1 + , 2 + emanating from some 
combination of the states of orbital symmetry [3 3 0] and 4 11] (these orbital states are 
degenerate in energy at B = 0).For B sufficiently negative, the nature of the ground state 
will again change. 

But let us focus on the region of B for which the triplet above is the ground state 
(B < —0.84). Besides the striking fact that the ground state is a triplet, what other 
evidence do we have of a 'phase transition' relative to the case where the orbital symmetry 
was [4 2 0]? Let us consider the case B = —1.0 and look again at tables I and 77 below the 
horizontal double lines. 

From table 777 we see that the main isovector E2 strength is from the J = + to the 
J = 2 + member of the ground state triplet i.e. a zero-energy transition. Although the 
B(E2) is substantial, —6.80 e 2 /m 4 , the rate would be zero if it is indeed a zero-energy 
transition. The isoscalar E2 strength is now split almost evenly between a low-energy state 
at 0.97 MeV and a high-energy state at 8.7 MeV. This is quite different from the case 
B > —0.67, where all the strength was concentrated in one state. 

Most interestingly a rather large isoscalar orbital Ml strength emerges (0.119 It is 
again a 'zero energy' transition, however from L = lS'=lJ = + toL=lS'=lJ = 1 + . 
Recall that for B = the isoscalar orbital strength is zero because the ground state has 
L = 0. The isovector scissors mode strength from J = + T = 1 to the J = 1 + T = 1 states is 
now fragmented into two parts -a low-energy part at 0.97 MeV with B(Ml) or i > i ta i = 0.13 fi 2 N 
and a high-energy part at 13.3 MeV with B(Ml) orUtal = 0.27 (i 2 N . For B > -0.67 all the 
1 + , T = 1 strength went only to one state. 

IV. CLOSING REMARKS 

We have shown that an isovector quadrupole-quadrupole interaction — xBQ(1)-Q(2)t(1)- 
t(2) with a large negative B yields very undesirable properties in p shell model calculations 
of 10 Be. On the other hand, such interactions appear to be needed to give correct splittings 
of the isovector and isoscalar giant quadrupole resonances. From our analysis of 10 Be in a p 
shell calculation, it appears that the best fit is obtained with B = (Elliott SU (3) model) 
-even a small positive B might be acceptable. 
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For B ~ —0.7 we get two sets of states collapsing to zero energy -first an isovector J = 2 + 
state, and then an L = 1 S = 1 triplet J = + , 1 + , 2 + . This behaviour is undesirable - 
no known nuclei behave in this way. Perhaps a remedy to this dilemma is to introduce 



momentum-dependent quadrupole terms in the interaction as Elliott had done Jll]. This 
enabled him to have an interaction which did not connect the AN = with the AiV = 2 
space. We would choose these so that the AiV = isovector quadrupole interaction is much 
weaker than the AiV = 2 part. 

Another suggestion is to bring effective mass into the picture when analyzing the sepa- 
ration of the isoscalar and isovector giant quadrupole resonances. The unperturbed energy, 
rather than being 2hu, is now 2^r. The energy of the isoscalar quadrupole resonance, rather 



than being \/2hu) is now \/2 hi f i IS|18|-f20|]. Using the 'empirical' energies and for 



(2C15 

v m * 



A3 A3 



the isoscalar and isovector quadrupole resonances respectively, we now modify the estimate 
of B in Sec. 2 as 

141 



63 -3 



With 2^ = 0.8 we get B = -0.97; with ^ = 0.7 we get B = -0.41. These are much 
smaller in magnitude than the value B = —3.1 we obtained with — = 1. Furthermore, in a 
AN = space, the magnitude of B will be even smaller due to renormalization effects ||. 
This argument is admittedly somewhat hybrid, but we believe it corresponds more closely 
to what happens when realistic interactions are used. 
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TABLES 





The energies 


ana n\hil, 


e p ,e n ) \ (in 


e 2 /m 4 ) to the first two 2 + states as 


a function 


of B 
















2j^ (isovector) 


a 


2^ (isoscalar) 


a 


B 


E (MeV) 




B(E2, 1, -1) 


E (MeV) 


B(E2, 1, 1) 


0.6 


3.49 




42.04 


4.02 


57.85 


0.4 


3.41 




39.02 


3.03 


63.75 


0.2 


3.00 




35.03 


2.48 


67.33 


0.0 


2.32 




33.40 


2.32 


65.09 


-0.2 


1.57 




30.26 


2.44 


67.71 


-0.4 


0.87 




29.49 


2.77 


66.45 


-0.6 


0.21 




29.11 


3.21 


65.06 



a The value of B{E2, 1, 1) to the state |2J) is zero. The value of B(E2, 1, -1) to the state \2%) is 



zero. 

TABLE II. The magnetic dipole orbital strengths (in jj? N ) -both isoscalar and isovector- as a 



function of B 







ISOSCALAR 




ISOVECTOR 








E 


1+, T = 1 


E 


1+, T = l(scissors) 


E 


1+, T = 2 




(MeV) 


B(M1) barbital 


(MeV) 


B(M1) Urbital 


(MeV) 


B(M1) Urbital 


0.6 


6.04 





9.27 


0.0040 


11.74 


0.0390 


0.4 


5.55 





8.59 


0.0043 


10.23 


0.0728 


0.2 


4.89 





8.34 


0.0370 


9.16 


0.1012 


0.0 


3.86 





8.49 


0.0895 (3^) 


8.49 


0.1492 ( j*) 


-0.2 


2.72 





8.93 


0.1431 


8.10 


0.1776 


-0.4 


1.64 





9.56 


0.1922 


7.91 


0.2019 


-0.6 


0.62 





10.31 


0.2314 


7.84 


0.2284 



TABLE III. Beyond the Cross Over Region, B = -1.0 



J = 2+, T = 1 

E*(MeV) 


B(E2, 1, 1) (e 2 fm 4 ) 


B(£2,l,-1) (e 2 /™ 4 ) 


0.00 





6.7120 


0.96 


15.93 





8.68 


21.68 





J = 1+, T = 1 
E*(MeV) 
0.00 
0.96 
13.3 


B(Ml) orbital ISOSCALAR (^) 
0.1194 




B(Ml) orbltal ISOVECTOR 


0.1368 
0.2672 
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Figure Captions 



Figure (1): The excitation energies of selected states in 10 Be as a function of the 
isovector quadrupole interaction parameter B. The solid line is for the 2\ state (isovector 
2 + ), the dashed line is for the 2 J state (isoscalar 2 + ) and the dot-dashed for the S = 1 L = 1 
triplet (J = 0+ 1+ 2+). 

Figure (2): Same as Figure 1 but for the J = 1 + T = 1 scissors mode branch (solid 
line) and for the J = 1 + T = 2 scissors mode branch (dashed line). 
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